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An ab initio study of magnetic exchange interactions in antiferromagnetic and strongly correlated 
3d transition metal monoxides is presented. Their electronic structure is calculated using the local 
self-interaction correction approach, implemented within the Korringa-Kohn-Rostoker band struc- 
ture method, which is based on multiple scattering theory. The Heisenberg exchange constants are 
evaluated with the magnetic force theorem. Based on these the corresponding Neel temperatures 
Tjv and spin wave dispersions are calculated. The Neel temperatures are obtained using mean field 
approximation, random phase approximation and Monte Carlo simulations. The pressure depen- 
dence of Tjv is investigated using exchange constants calculated for different lattice constants. All 
the calculated results are compared to experimental data. 



I. INTRODUCTION 

In the last years there has been a general strong inter- 
est in finding materials with specific or even parametris- 
ablc magnetic properties. Such materials could be useful 
in the field of spintronics. A lot of the promising can- 
didates are strongly correlated electronic systems which 
in many ways are still a challenge to be properly de- 
scribed theoretically regarding their electronic ground 
state properties. On the other hand, for reliable predic- 
tions about magnetic properties of materials it is essential 
to have theories describing the magnetism adequately by 
quantitative and qualitative means. One of these theories 
is the Heisenberg theory of magnetism, which we shall ap- 
ply in the present paper. Its central quantities are the 
Heisenberg exchange constants Jy , which are of general 
fundamental interest. In particular they provide infor- 
mation about the magnetic periodicity (via their Fourier 
transform), the spin- wave dispersion, magnetic critical 
temperatures and also allow predictions on structural ef- 
fects caused by magnetism 1 ^. 

In this paper we concentrate on the study of the mag- 
netic exchange interactions of transition metal monoxides 
(TMOs), specifically MnO, FeO, CoO and NiO. They are 
charge-transfer insulators, well known for strong correla- 
tion effects associated with the TM 3<i-electrons. Orig- 
inating from the Anderson-type superexchange, their 
cqulibrium magnetic structures are of the antiferromag- 
netic II (AFII) order, characterized by planes of opposite 
magnetization which are stacked in (11 Indirection. Re- 
cently, a Mott transition has been observed in MnO at 
high pressure of about 105 GPa, in resistivity 3 - and X- 
ray spectroscopy measurements* 4 ^ which stimulated new 
theoretical studies in this high pressure region^ There 
already exists a large body of neutron scattering measure- 
ments of magnetic structures and magnetic excitations in 
transition metal monoxides. However, the development 



of new experimental techniques such as neutron pow- 
der diffraction^uiiiii and polarized neutron reflectivity^ 1 - 
has renewed interest in studying TMOs as antiferromag- 
netic benchmark materials. Modern neutron spectrom- 
eters operate with such a high efficiency that also high 
angle diffraction experiments can be performed to un- 
ravel complex magnetic order e. g. in thin films 1 - 2 -. 
From the theory point of view conventional methods 
such as the local spin density approximation (LSDA) to 
density functional theory (DFT), treating electron cor- 
relations at the level of the homogeneous electron gas, 
fail to provide an adequate description of the electronic 
structure of these oxides. Over the years a number 
of approaches has been developed, aiming at improve- 
ments to the LSDA treatment of electron correlations, 
and applied to TMOs with varying degrees of success. 
Among them arc: the LSDA+U method 1 ^ 4 -, GGA+Ui^, 
self-interaction corrected (SICj-LSDAi^^i!^^, hy- 
brid functional a 22 i 23 , and finally dynamical mean field 
theory^ 4 -. In general, they have improved lattice con- 
stants, band gaps and magnetic properties, some of them 
have also obtained good agreement with spectroscopics. 
In the present paper we shall use the so-called local self- 
interaction correction 2 ^ (LSIC) scheme for the calcula- 
tion of the electronic ground states of the TMOs. As the 
aim is the investigation of magnetic interactions we com- 
bine the LSIC scheme with the magnetic force theorem 
(MFT)2£ in order to obtain the Heisenberg exchange pa- 
rameters Jij. 

The LSIC scheme 2 ^ is based on the implementation of the 
SIC-LSDA formalis m 16 i 17 i 25 within multiple scattering 
theory in the framework of the Korringa-Kohn-Rostoker 
(KKR) band structure method. It was first applied to 
/- electron system o 27 i 28 , but recently also to TMOs2i. 
Within the KKR method one calculates the Green func- 
tion of the investigated system. This Green function 
is then straightforwardly used in the application of the 
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MFT. It is for the first time that this combined approach 
is applied for calculating exchange constants of the tran- 
sition metal monoxides. The results of that are compared 
to the exchange constants extracted from the total en- 
ergy differences for a number of magnetic structures and 
mapping them onto a Heisenberg Hamiltonian. Most of 
the earlier applications of the latter approach have been 
based on the assumption that only the first two exchange 
interaction constants arc nonzero. Although the present 
combined approach also relies on the mapping onto a 
Heisenberg Hamiltonian the assumptions regarding the 
number of non-zero exchange constants are not needed, 
which is advantageous to systems with reduced symme- 
try such as thin films and layered structures (where the 
justification for such an assumption is not clear from the 
very beginning). 

Having calculated the for the ground states of the 
TMOs we also calculate and discuss them as a function 
of external pressure for moderate values of the latter. 
This is mainly inspired by the recent high pressure mea- 
surements of TMOs^£. 

Based on the calculated magnetic exchange interactions 
the transition temperatures can be obtained. In this pa- 
per it is done in three different ways, namely by applying 
mean- field approximation (MFA) , random phase approx- 
imation (RPA), and using classical Monte Carlo (MC) 
simulations. The respective results are then compared to 
those obtained from the disordered local moments (DLM) 
method^, which does not involve mapping onto a Heisen- 
berg Hamiltonian, but is based on the same ground state 
electronic structure calculations as the present paper—. 
The last subject we focus on are magnetic excitations. 
On one hand, with given Jy , one can calculate the 
magnon spectrum of any material. On the other hand, 
measuring the latter experimentally is a direct method 
to examine its exchange constants. Thus, comparing cal- 
culated to experimental spin wave dispersions provides a 
straightforward tool for determining the accuracy of the 
calculated Jy. 

The present paper is organized as follows: In section 
HT1 the theoretical approaches for the calculation of elec- 
tronic structure, exchange interactions and Neel temper- 
atures are presented. The computational details are de- 
scribed in section HTIl Scction HVl contains the results and 
discussion. The exchange parameters and the Neel tem- 
peratures are presented for theoretical equilibrium lattice 
constants and as a function of lattice constants and pres- 
sure, respectively. Finally, the calculated magnon spectra 
of the TMOs are discussed in reference to experiments. 
The paper is concluded in section fVl 



II. THEORY 
A. Electronic Structure 

For the electronic structure calculations of TMOs we 
use a multiple scattering theory-based implementation of 



the SIC-LSDA method^ 5 ., whose total energy functional 
is 

^SIC-LSDA^^ = E LSBA [nhni] _ 

^(£ H [n a ] + £ x L c SDA [^-0]) ,(1) 

with the LSDA energy functional in units of Rydberg 
given by 

E LSDA [n hni ] = ]T <</w|-V 2 \<j> aa ) + 

E cxt +E H [n}+E^ DA [n h n L ] . (2) 

Here (j) aa is a Kohn-Sham orbital, aa a multi-index la- 
belling the orbitals and spin (t or J,), respectively, n ac = 
\4>aa\ 2 , n a = X)a n acr and n = Tif + 7i|. £; LSDA differs 
from E LSI)A since the kinetic energy is evaluated with re- 
spect to the orbitals minimizing the SIC-functional. The 
summations run over all the occupied orbitals, E cxt de- 
notes the external energy functional due to ion cores, £"h 
is the Hartrce energy and i4 SDA is the LSDA exchange- 
correlation energy functional. The second term in Eq. (fT]) 
is the so-called self-interaction correction^ 5 , for all the oc- 
cupied orbitals a. It restores the property 

E B [n a ]+E^ act [n a . a ,0} = 0, (3) 

that the exact DFT exchange-correlation functional has, 
namely that for any single orbital density the Hartree 
term should be cancelled by the corresponding exchange- 
correlation term. The cost paid for restoring the above 
property is the orbital dependence of the SIC-LSDA en- 
ergy functional (Eq. (TT|)). The correction is only substan- 
tial for localized orbital states, but vanishes for itinerant 
states. In the limit of all itinerant states the SIC-LSDA 
total energy functional is identically equal to the LSDA 
functional. 

The main idea behind the "local" implementation of the 
SIC-LSDA formalism (LSIC) is that within multiple scat- 
tering theory, in the framework of the KKR method, one 
works with the scattering phase shifts, describing scat- 
tering properties of single atoms in a solid. Among them 
only the resonant phase shifts are relevant, as they refer 
to localized states. Thus the self-interaction correction is 
associated with the on-site scattering potentials and leads 
to modified resonant scattering phase shifts. In particu- 
lar, they become stronger localized. Details of the LSIC 
implementation are discussed in Ref. l27l 

B. Magnetic Interactions 

The Heisenberg theory of magnetism assumes that it 
is possible to map magnetic interactions in a material 
onto localized spin moments, which in a classical picture 
can be represented by a vector. The resulting classical 
Hamiltonian, 

H = -~y] J l] e l ■ ej , (4) 
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contains only the unit vectors e^-j of the spin moments 
and the exchange parameters describing the interac- 
tions between them^i. Here i and j index the sites. 
It should be mentioned here that the Hamiltonian (fJJ 
can be extended to include additional effects like magne- 
tocrystalline anisotropy or tetragonal or rhombohedral 
distortions of the lattice. The latter reflect magnetoe- 
lastic effects which result in two different values for the 
nearest neighbour (NN) exchange parameters, depending 
on a parallel or antiparallel alignment of the moments. 
Such effects are usually present in experiments. Thus, 
special care is required when comparing theoretical and 
experimental results. 

Our method of choice for the calculation of the exchange 
parameters Jij makes use of the magnetic force theorem, 
but invokes also mapping onto a Heiscnbcrg Hamiltonian. 
For comparison, we also apply the most commonly used 
approach which relies on the calculation of total energy 
differences between different magnetic configurations and 
mapping them onto a classical Heisenberg Hamiltonian. 

1. Magnetic Force Theorem Approach 

The idea behind the magnetic force theorem^ is to 
consider infinitesimally small rotations of classical spins 
at two different lattice sites. These give rise to energy 
changes that are mapped onto the classical Heisenberg 
Hamiltonian via multiple scattering theory. This ap- 
proach is based on the assumption that the potentials 
arc unchanged by the rotations. The advantage of the 
MFT method is that the exchange integrals can be cal- 
culated directly in the relevant magnetic structure. The 
result for the exchange parameter of the two magnetic 
moments at sites i and j can be written as 

J y = ^ / " deImTr L (A 4 ff Aj-rf + A^A.ff ) , 

( 5 ) 

where f u is the scattering path operator between sites i 
and j and A, = — ^Ij. 1 ^ with h being a single scattering 
operator for the atom at site i. If not stated otherwise, 
all the results discussed later would have been obtained 
with the exchange parameters calculated using Eq. ([5]). 

2. Energy Differences Approach 

In this approach the total energies of the TMOs in the 
ferromagnetic (FM) and antifcrromagnctic I and II (AFI 
and AFII) configurations are taken into account. The 
AFI structure is characterized by oppositely magnetized 
planes which are stacked in (lOO)-direction. Suppose that 
magnetic interactions operate only between TM atoms 
- an assumption which is to be discussed later — the 
mapping onto the Heisenberg Hamiltonian yields 

Ji = 4 (^afi - E FM ) (6) 



and 

^ 2 = 7Z (4^AFii — 3E AF i — E FM ) , (7) 

where J\ describes the interaction between the NN and 
J2 that of the next nearest neighbours (see Fig. [TJ. This 
mapping also assumes that the interaction between NN is 
independent of the sublattice the TM atoms are located 
on. Of course the choice of the three above mentioned 
structures restricts one to the determination of Jj and 
J2 only. Using more magnetic structures and hence cal- 
culating more exchange parameters is in principle possi- 
ble. However, due to the nature of the present exchange 
mechanism, the super exchange, this has usually not been 
done for the TMOs. Although with this method we also 
restrict exclusively to J\ and J2 , it is hoped that the com- 
parison with the MFT method will shed some light on the 
validity of the underlying assumptions for the TMOs. 




FIG. 1: (color online) Schematic representation of the mag- 
netic interactions in a (100) plane of the rocksalt structure of 
TMOs. The TM ions (blue) interact via Ji (dotted arrows) 
with their nearest and via J2 (solid arrow) with their next 
nearest neighbours. In Anderson's super-exchange picture the 
indirect exchange is mediated by O ions (yellow circles), re- 
sulting in a 90°- and a 180°- exchange interaction for Ji and 
J2, respectively. Note that Ji also contains contributions from 
direct overlap. 



C. Neel Temperatures 

Having calculated the exchange parameters for an an- 
tiferromagnet one is able to calculate the Neel tempera- 
tures, Tn. Several different approximations can be used, 
in particular the mean field approximation, the random 
phase approximation and classical Monte Carlo simula- 
tions. 
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1. Random Phase Approximation 

In the RPA one solves the equation of motion for the 
Green function of the spin operators. Following the ap- 
proach of Rusz et a/i^i, one ends up with a semi-classical 
formula for the average spin polarization (e z A ) (T) of a 
sublattice (e.g. A), function of temperature T, 



(el) (T) = £ 



~|dq [N^q)] 



,1,1 



(8) 



Here C(x) is the Langevin function, C(x) = coth(a;)— 1/x, 
is the volume of the first Brillouin zone and kg denotes 
the Boltzmann constant. The matrix elements of N(q) 
are defined as 

Nab(ci) = S AB Jac(0) (e z c ) ~ (e z A ) Jab(<1), (9) 



with the Fourier transforms of the exchange parameters 
given by 

Jab (q) = IJ2 Jike^-^hflu , (10) 

where v denotes the number of interacting magnetic sites 
and 7^ equals one if site j is on the magnetic sublat- 
tice A and zero otherwise. Eq. (|5|) has to be solved self- 
consistcntly since the unknown quantity appears on its 
left and implicitly also on the right, via N(q). The Neel 
temperature is equal to the highest value of T at which 
(e A ) (T) becomes different from zero. 



2. Mean-Field Approximation 

To obtain the MFA estimate of the Neel temperature, 
a matrix with elements 



Qab = ^-Jab(0) 



(11) 



is constructe d 32 ' 33 , where Jab(0) stands for the Fourier 
transform of the exchange parameters, defined via 
Eq. (fTU|) . at q = 0. The largest eigenvalue of yields the 
Neel temperature. If for any TMO in the AFII structure 
only the nearest and next-nearest neighbour interactions 
are considered (Ji and J2, respectively), then the largest 
eigenvalue yields the well-known relation, /sbTn = 4J2, 
indicating that the nearest neighbour interaction J\ does 
not have any influence on T^. Since fluctuations are com- 
pletely neglected in MFA the resulting Neel temperatures 
are commonly overestimated. 



3. Monte Carlo Simulations 

We give a rather brief summary of the method of MC 
simulations as they arc performed in this paper. For a 



deeper and complete understanding we refer the reader 
to the book by Landau and Binder—. 
To estimate T/v via MC simulations a lattice representing 
the structure of the investigated system is constructed. 
The magnetic moment at lattice site i interacts with its 
neighbours j via the Jij . During a MC run one picks a 
lattice site j with the magnetic moment vector Bj, cre- 
ates a new random direction e^- and decides by looking 
at the energy of the system whether e^- is accepted or 
is kept. Performing this procedure A-times on a lattice 
of N sites is defined as one MC step. 
Starting from a certain initial configuration the system is 
brought into thermal equilibrium for a fixed temperature. 
After this, "measurements" and thermodynamical aver- 
aging of the observables of interest are performed. Since 
it is impossible during a simulation run to go through all 
possible configurations of the system, which would be for- 
mally necessary for averaging, one must ensure that the 
configurational subspace that one is restricted to is of 
physical significance. This is done by performing the so- 
called importance sampling. It is applied when one has 
to decide between the old and new magnetic moment 
vectors e'j and Bj described above. There exist several 
methods to do this, in the present paper the Metropolis 
algorithm 3 ^ is used. 

One must be aware of the fact that the finite size of the 
lattice, despite being periodic in all 3 dimensions, leads 
to a systematic error in the determination of the crit- 
ical temperature. This so-called finite-size effect, how- 
ever, becomes smaller with increasing lattice size. It can 
therefore be eliminated by extrapolating the critical tem- 
peratures for different lattice sizes. 

For a magnetic system as in the present case it is straight- 
forward to measure two quantities. One is the staggered 
magnetization^ m s being some sort of an average of the 
absolute values of magnetization of the sublattices, 



1 N 

-Y 



(12) 



Here, j labels the lattice sites being N in total, Bj is the 
unit vector of the magnetic moment at lattice site j and 
Q is the normal vector of the planes of equal magneti- 
zation, in the AFII structure being (1, 1, 1) for example. 
As a note, m s is used instead of the total magnetiza- 
tion since the latter is equal to zero in antiferromagnets. 
The other quantity measured is the inner, i.e. magnetic, 
energy of the system E, which is given by Eq. Q. This 
whole procedure of relaxing into thermal equilibrium and 
thermodynamical averaging is repeated for different tem- 
peratures. In principle one can determine T/v from the 
slope of the temperature dependence of m s (T) and E(T). 
However, there are quantities that show the critical tem- 
peratures more clearly. These are the magnetic suscep- 
tibility derived from m s , the specific heat derived from 
E, and the 4th-order cumulam— U4. The first two have 
a singularity at T = Tjy, the last one has the property 
that the curves Kt(T) calculated for different lattice sizes 



5 



have a crossing point at T = TV- 

4- Quantum Effects 

As already mentioned, the three approaches described 
above are based on mapping of the single magnetic mo- 
ment interactions onto a classical Heisenberg Hamilto- 
nian. These moments, however, are quantum objects and 
this should in some way be accounted for in the calcula- 
tions. Wan, Yin and Savrasov^ and Harrison 3 ^ did this 
by replacing the classical S 2 in the Heisenberg Hamil- 
tonian with the quantum mechanical expectation value 
S(S + 1), when calculating magnetic properties. Since 
in Eq. (J3]) S 2 is included in the Jij— , then to be consis- 
tent, one has to divide again by S 2 . This gives rise to 
a factor (S + l)/S for the energy and, eventually, also 
for the Neel temperature 4 ^ This factor is, however, not 
needed when the Neel temperature is obtained based on 
the DLM method since it does not explicitely use the Jij . 

D. Magnon Spectra 

With the exchange interactions determined one can 
also calculate the magnon spectra E(q). They are of 
special interest since they provide the standard method 
for determining the exchange parameters experimentally. 
The latter would be done by fitting a Hamiltonian con- 
taining the as fitting parameters to a measured spin 
wave dispersion. As already mentioned in section[lll such 
Hamiltonians usually contain more terms than the one 
given in Eq. (j4]), which is why comparisons between dif- 
ferent Jij results must be done carefully. 
Considering multiple sublattices one can define magnon 
spectra as the eigenvalues of the matrix N(q) given 
by Eq. ([9]). Assuming two magnetic sublattices, with 
the same absolute magnetization, and considering only 
the nearest and next-nearest neighbour interactions, the 
spectra are given b y 30 ' 41 

E(q) = - ^(J ++ ( q )-H Q f-J 2 _(q) . (13) 

Here \i is the magnetic moment of the two sublattices in 
units of /IB, J++(q) (J+_(q)) are the Fourier transforms 
of the intra- (inter-)sublattice exchange parameters ex- 
pressed respectively as 

J ++ (c l )=2J 1 X (14) 
x (cos7ra(q y + q z ) + cosira(q x + q y ) + cos-n:a{q x + q z )) 

and 

J+_(q) = 2J lX (15) 
x (cos7ra(g. y - q z ) + cosira(q y - q x ) + cosira(q x - q z )) 
+ 2J2 (cos 2na q x + cos 2ira q y + cos 27ra q z ) , 

and H = J ++ (0) - J+_(0) = -6J 2 . In Eqs. {HJ) and 
(I15p the vector q and accordingly its components are 



needed in units of 2n/a with a being the lattice constant 
of the TMO considered. 



III. COMPUTATIONAL DETAILS 

The transition metal monoxides crystallize in the rock- 
salt structure (Bl, Fm3m, space group 225). At low 
temperatures they show small lattice distortions (< 2%). 
However, these distortions are not considered in the 
present calculations. The crystal potentials for the 
ground state calculations are constructed in the atomic 
sphere approximation (ASA) . The ASA radii for the TM 
and oxygen atoms are chosen as 0.2895 a, with a being 
the lattice constant of a given TMO. To reduce the ASA 
overlap while keeping a good space filling, empty spheres 
are used with the ASA radii equal to 0.1774 a. The ratios 
of the respective ASA radii are kept constant across the 
TMO series. 

For the electronic structure calculations the complex en- 
ergy contour has 24 Gaussian quadrature points, and for 
the Brillouin zone (BZ) integrations a 14x14x14 k-points 
mesh is constructed. For the calculation of the mag- 
netic interactions, using the MFT, 60 energy points on a 
Gaussian mesh in the complex plane are chosen. Conver- 
gence of the Jij with respect to the number of k-points is 
achieved with a 20x20x20 k-points mesh per energy point 
for the first 50 of them, and a 60x60x60 k-mesh for the 
last 10 energy points, lying close to the Fermi energy. 
For the MC simulations an fcc-lattice representing the 
transition metal atoms in the TMO crystal is con- 
structed. To avoid finite-size effects, the size of the lat- 
tice is varied from 40x40x40 to 60x60x60 elementary fee 
cells. To use all observables described in the MC part 
of Section |TT] one has to restrict the simulations to a rel- 
atively small number of MC steps. This is necessary in 
order to prevent the system from changing the orienta- 
tion of the ferromagnetic sublattices, for example from 
(111) to (111), which are degenerate in energy. Thus, 
starting from the AFII state, the system is assumed to 
have reached thermal equilibrium after 5,000 MC steps, 
and for averaging 10,000 MC steps are performed. If 
one does so all observables yield the same result for the 
Neel temperatures for each TMO, respectively. To en- 
sure a thorough exploration of phase space, simulations 
with up to 100,000 MC steps for averaging have also been 
performed. In this case the specific heat, not affected by 
reorientations of the magnetic sublattices, has reassur- 
ingly indicated the magnetic phase transitions to occur 
at the same temperatures as in the short simulations. 



IV. RESULTS AND DISCUSSION 

A. Exchange parameters 

The present calculations of the exchange parameters of 
TMOs use the ground state electronic structure proper- 
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ties of these materials as input. The latter are obtained 
self-consistently with the LSIC method, explained in de- 
tail in Ref. Hi]. In particular, as seen in Eq. ((5]), for the 
MFT approach the relevant quantities are the scattering 
properties evaluated at the equilibrium lattice constant 
of the ground state, AFII, magnetic structure. For the 
approach based on the energy differences only the to- 
tal energies of the FM. AFI and AFII structures, eval- 
uated at the theoretical equilibrium lattice constants of 
the AFII configuration, are of relevance. 

TABLE I: The calculated (calc.) equilibrium lattice constants 
and spin magnetic moments, per TM-atom, for all the stud- 
ied TMOs in AFII structure. The oxygen atoms are not spin 
polarized in the AFII environment and the induced moments 
on the empty spheres are very small. Consequently the calcu- 
lated spin magnetic moments of TMOs are practically equal 
to those of their TM-atoms. The experimental (exp.) values 
of the magnetic moments contain not only the spin but also 
orbital contribution. 



TMO 


do 


[A] 


/i [hb] 


calc. 


exp. 


calc. exp. 


MnO 


4.49 


4.44« 


4.63 4.54^ 


FeO 


4.39 


4.33^. 


3.68 3.32^ 


CoO 


4.31 


4.26^ 


2.69 2.40^ 


NiO 


4.24 


4.17a 


1.68 1.90^ 



From Table U we can sec that the LSIC method, 
treating localized and itinerant electrons on equal 
footing, reproduces well the equilibrium lattice constants 
and also the corresponding spin magnetic moments in 
the AFII structure. The overall agreement with the 
experimental values is reasonable for both quantities. 
Note, however, that the experimental magnetic moments 
listed in the table include both the spin- and orbital- 
contributions, which are substantial for FeO and CoO, 
and non-negligible even for NiO. Regarding the calcu- 
lated spin magnetic moments, they are effectively equal 
to the spin moments of the TM atoms, as the oxygen 
atoms are not polarized in the AFII environment, and 
the induced spin moments on the empty spheres are very 
small. In addition, as seen in Fig. [21 the spin magnetic 
moments show considerable dependence on the lattice 
constants, indicating that a similar behaviour may also 
be expected for the calculated exchange constants. This 
dependence of the calculated spin magnetic moments 
on the lattice constant was also observed in previous 
studies ^ 15 ' 48 and its trend agrees with the Stoner model, 
stating that magnetic moments eventually collapse at 
very high pressures. 

Using the above ground state properties in the MFT 
approach, we have calculated the exchange constants 
for the first 11 neighbour shells. As expected only 
the first two of them, J\ and J2, are of relevance as 
those corresponding to the higher shells are less than 
0.1 me V in magnitude. This agrees well with the idea 
of super-exchang and can easily be explained with 
it. Consequently, and for the purpose of comparison 




a [A] 



FIG. 2: (color online) The dependence of the calculated spin 
magnetic moments on lattice constants. The spin magnetic 
moments have been divided by their respective equlibirium 
values given in Table [I] The crossing points of each curve with 
the horizontal dotted line at fij ' fj, eq — 1 mark the equilibrium 
lattice constant of each TMO. 



with the results of the energy difference approach, in 
Table |TT] we display only the J\ and J2 quantities, 
as well as the experimental results. Our results also 
provide justification for one of the assumptions underly- 
ing Eqs. ((6]) and ((7|), that the interactions between NN 
atoms are not dependent on the sublattices of the atoms. 

From our results in Tabic ITT1 one finds that the J2 pa- 
rameters constitute the major part of magnetic exchange 
in TMOs and that in magnitude they agree reasonably 
well between the two theoretical approaches, MFT and 
AE. The results are about 70 to 80 % of the experimen- 
tal values, except for FeO, where the agreement for the 
MFT J2 is almost perfect and the AE value is larger than 
the experimental one. This rather accidental agreement 
for FeO can most likely be attributed to the fact that the 
experimental values are measured for Wiistite samples 
Fei- x O with x ^ 0.— Also, the experimentally observed 
trend of the increasing absolute value of J2 across the se- 
ries is present in both approaches and is most likely asso- 
ciated with the increasing number of the TM 3d electrons, 
responsible for the magnetic super-exchange. Regarding 
the quantitative agreement with experiment, the MFT 
results agree better on average. One could envisage that 
this agreement could be further improved, if the MFT ap- 
proach was applied in the DLM stat o 29 ^ 56 . Nevertheless, 
compared to the other calculations displayed in Tables [TT1 
and IIIII our present results may already be considered as 
being at least as good as those. 

The situation is very different for the J\ exchange param- 
eters. From Table [TT1 we see that, with the exception of 
NiO, the absolute magnitudes of J\ are, within about 
30%, similar between the two theoretical approaches, 
but the signs are opposite. Looking more closely at 
the parameters calculated with the MFT we see that 
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TABLE II: The exchange parameters, J;, in meV, with i being the shell index, for the first two shells of all the studied TMOs 
and both the MFT and AE, energy difference, methods. The experimental (exp.) values are given in the leftmost column, 
respectively, according to the Hamiltoniar*^ in Eq. Q. For i > 2, the absolute magnitudes of the J;'s for all the TMOs have 
been less than 0.1 meV. For both the MFT and AE approaches the calculated equilibrium lattice constants from Table|T]have 
been used. The results of Harrison^, obtained using a tight-binding (TB) formalism for the complete series of TMOs are listed 
here for a direct comparison. The other previous results, obtained for selected monoxides, are listed in Table ITlTI 





exp. 


MFT 


AE 


TB 39 


exp. 


MFT 


AE 


TB 39 


MnO 


-2.06,-2.64^ 


-0.91 


0.68 


-4.41 


-2.79^ 


-1.99 


-1.65 


-1.09 


FeO 


1.04, 1.8422, 


-0.65 


0.48 


-2.99 


-3.24^ 


-3.17 


-3.50 


-1.56 


CoO 


0.70^-, -1.07-2 2 - 


-0.32 


0.53 


-1.83 


-6.30£i, -5.31-22 


-4.84 


-4.40 


-1.64 


NiO 


-0.69^,0.69^ 


0.15 


1.42 


-1.44 


-8.66^,-9.51^ 


-6.92 


-6.95 


-1.88 



TABLE III: Summary of the first principles results for Ji and J2 in MnO and NiO, based on the Hamiltonian^i in Eq. @, 
from the present and previous theoretical works for comparison. Only those close to experimental values are listed. For details 
see the corresponding references. We found one result by Feng2£ for CoO, obtained by using the B3LYP hybrid functionals 
method, Ji = —47.12 meV and J2 = —42.56 meV. To our knowledge no further theoretical papers giving numerical values for 
the exchange parameters of FeO exist. 





MnO 






NiO 




method 


Ji [meV] 


J2 [meV] 


method 


Ji [meV] 


J 2 [meV] 


exp. 


-2.06,-2.64 


-2.79 


exp. 


-0.69,0.69 


-8.66,-9.51 


this work 


-0.91 


-1.99 


this work 


0.15 


-6.92 


LDA+U^I 


-2.50 


-6.60 


GGA+U±2 


0.87 


-9.54 


OEPii 


-2.85 


-5.50 


SIC-LMT02£ 


0.90 


-5.50 


PBE+U^ 


-2.21 


-1.16 


Fock35^ 


0.95 


-9.35 


PBEC 22 


-3.10 


-3.69 


B3LYR22 


1.20 


-13.35 


HF-22. 


-0.73 


-1.16 


UHF— 


0.40 


-2.30 


B3LYR22 


-2.64 


-5.52 









the results show the opposite trend to that found for 
Ji. Namely, the antiferromagnetic coupling is getting 
weaker as one moves from MnO to CoO, and in NiO it 
becomes ferromagnetic. This can be explained by as- 
suming both kinds of interaction to be present and to be 
competing, in the direct and indirect exchange between 
NN TM atoms. The picture is relatively intuitive for 
the direct case. For Mn, which has half filled d-shclls, 
one expects antiferromagnetic coupling since an electron 
hopping from one Mn atom to the other one keeps its 
spin. Thus, this transfer clearly prefers antiferromag- 
netic alignment of the Mn atomsS. Moving across the 
TMO series the occupation of the minority spin channels 
is growing. This increases the probability of an electron 
hopping, if the TM atoms are ferromagnetically aligned. 
Thus, the character of the exchange should go towards 
ferromagnetic, which is what we find for the J\ calcu- 
lated via the MFT. Regarding the indirect exchange, we 
can also follow GoodcnougbJii. Nearest TM neighbours 
interact antiferromagnetically when two electrons in the 
same oxygen p- or s-orbitals are excited to the empty TM 
e g -orbitals. The strength of this kind of interaction can 
be assumed not to change a lot along the TMO series, 
since the occupation of both the oxygen p- or s- and the 
TM e g -orbitals does not change either. Ferromagnetic 



coupling on the other hand is provided by electrons of 
alike spin that are in different orbitals of the O atom, ft 
is strengthened by a growing occupation of i2 g -orbitals 
because this increases intraatomic exchange. Since the 
fog-occupancy is rising when moving across the TMO se- 
ries from MnO to NiO^i one would expect that the mag- 
nitude of the ferromagnetic interaction increases while 
the antiferromagnetic does not. This tendency is clearly 
present in the MFT values for J\ in Table HT1 
Looking at the agreement with experiment not much 
overlap can be spotted. For the MnO the agreement 
is satisfactory considering the simplicity of our Hamilto- 
nian. For FeO the sign is opposite. This could be caused 
by the above mentioned fact that in experiment Wiistite 
samples of the kind Fei-^O are investigated while our 
calculations are performed for the ideal FeO system. For 
CoO and NiO comparison is difficult since experimental 
values of opposite signs, but similar absolute magnitude 
have been measured in different experiments. To con- 
clude the comparison of MFT-Ji and experimental J\ 
one can say that the agreement is not as good as for the 
J 2- Possible reasons have just been given, but it also 
seems that the experimental dctcrminination of the J\ is 
not as accurate as for the Ji , as can be seen from the va- 
riety of numbers obtained for the same compounds. The 
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lack of agreement between experimental and MFT values 
for the Ji may influence the calculations based on those. 
For the magnon spectra it could be expected that, be- 
sides quantitative differences, due to the different signs 
even qualitative changes of the curves might occur. Wc 
shall see later, however, that the latter is not the case. 
The effect on the calculated Neel temperatures should 
be small in any case since the energy contributions of the 
NN in the AF11 structure are canceled out. 
For the J\ calculated with the energy difference approach 
no obvious trends are seen in Table UH and in addi- 
tion they are positive for all TMOs. For NiO the latter 
agrees qualitatively with the MFT- Ji and also with pre- 
vious theoretical results, and the agreement with those by 
Kodderitzsch et aZ.— , J\ = 0.9 meV and J2 = —5.5 meV, 
is also quantitatively rather good. For the other TMOs 
the sign of J\ is opposite to the ones calculated with the 
MFT, and for MnO and CoO they also do not agree with 
previous theoretical investigations. The totally different 
behaviour — compared to the MFT-values — can be ex- 
plained by looking at the electronic ground states of the 
calculated AFI and FM structures. In both of them the 
oxygen atoms carry a magnetic moment, which they do 
not in the AFII structure. This magnetic moment can 
be assumed to give rise to magnetic interaction with the 
neighbouring TM atoms (as a matter of fact, applying 
the MFT to AFI or FM structures yields NN exchange 
parameters of several meV in magnitude between TM 
and oxygen atoms). This, however, is in strong contra- 
diction to the assumptions underlying Eqs. ((6J and (|7|). 
stating that magnetic interaction only occurs between 
TM atoms. Thus, when using these equations anyway, 
this "artificially" created magnetic exchange is projected 
onto the J\ and J2. The reason for the latter quantity 
being relatively close to its counterpart calculated with 
the MFT is probably due to the large energy differences 
between the AFII configuration and the AFI as well as 
the FM configuration for each of the TMOs. This obvi- 
ously reduces the error made in Eq. ([7]). 
It should be mentioned that for all calculated pairs of J\ 
and J2, using SIC-LSDA, the resulting ground state mag- 
netic structure is that of AFII^, despite the relatively 
large spread of the J\ parameters. Note that the ex- 
change parameters J\ and J2, obtained for the TMO se- 
ries by applying the MFT approach to the LSDA ground 
state electronic structure, show no agreement with exper- 
iment, except for NiO, which can perhaps be considered 
as a lucky coincidence. Furthermore, the J's are longer- 
ranged, i.e. their character is more metallic. This agrees 
with the fact that their uncorrected (no SIC) ground 
states show only very small or no band gaps at all^I. 
Finally, we would like to comment on the variation of 
the exchange parameters as a function of lattice constants 
shown in Fig.[3]for all the TMOs. As one can see, the ab- 
solute value of J2 increases with decreasing a. This is in 
good agreement with the interpretation of the exchange 
parameters in terms of overlap integrals. The closer the 
atoms are, the larger the overlap is between the TM d- 



4.0 4.1 4.2 4.3 4.4 4.5 4.6 




4.0 4.1 4.2 4.3 4.4 4.5 4.6 
a [A] 



FIG. 3: (color online) (a) The Ji- and (b) ^-dependency on 
the lattice constant o for all TMOs, calculated using MFT. 
The vertical lines mark the calculated equilibrium lattice con- 
stants from Table U 

orbitals and the oxygen p-orbitals. A similar behaviour 
is found for the J\ , which can be understood in the same 
way as for J2. However, going through the TMO se- 
ries and starting with MnO, the change of the J\ gets 
smaller as the antiferromagnetic character becomes less 
pronounced. According to Goodenough's arguments^, 
this suggests that the ferromagnetic coupling becomes 
more prominent than the antiferromagnetic one. 



B. Neel temperatures 

The calculated transition temperatures are summa- 
rized in Table IIVI One finds that MFA overestimates 
the experimental Neel temperatures, whereas RPA un- 
derestimates them. This is what can be expected from 
general considerations^. One can also see that the Neel 
temperatures calculated in the RPA approach based only 
on Ji and J2 do not differ significantly from those cal- 
culated using the II neighbour shells. This again agrees 
with the idea of superexchange. What is not expected is 
that, for MnO and FeO the RPA and MC results are rel- 
atively small compared to experiment. In fact, for these 
TMOs the MFA gives a better estimate. The probable 
reason for that is the general relative underestimate for 
the Ji- The latter being the main contribution of mag- 
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TABLE IV: Summary of the Neel temperatures calculated 
with the Jij from the MFT approach (sec Tablc Ul) . In the top 
two rows the experimental and the DLM values are listed, fol- 
lowed by the RPA values based on the interaction of the first 
11 TM-TM-shells and of only the nearest and next-nearest 
neighbours (i.e. only Ji and J2). In rows 5 and 6 the MFA 
results shown, again using 11 or 2 shells, respectively. In the 
last row the results of the Monte Carlo simulations are pre- 
sented. 



T N [K] 


MnO 


FeO 


CoO 


NiO 


Experiment 


118 


192 


289 


523 


DLM22 


126 


172 


242 


336 


RPA with Ji_n 


81 


146 


252 


440 


RPA with Ji, 2 


87 


155 


260 


448 


MFA with Ji_n 


122 


210 


362 


628 


MFA with Ji, 2 


129 


221 


373 


644 


MC 


90 


162 


260 


458 




10 20 

P [GPa] 

FIG. 4: The normalized RPA-based Neel temperatures 
TV, RPA rpa f° r a U TMOs as a function of pressure p. 
Tn, rpa i s taken from row 3 of Table ITVl 



netic exchange, their underestimate is largest for MnO 
and decreases towards NiO. An exception is FeO. The 
agreement for J 2 is almost perfect, yet the RPA and MC 
estimates are roughly of the same quality as those for 
the other TMOs. However, it can again (see discussion 
of Ji and J2) be argued that due to the experimental 
imperfect FeO lattice other effects not considered in our 
approach may play an important role for the formation 
of magnetic order. The DLM results of Hughes et al2& 
are, with the exception of NiO, in good agreement with 
the experimental values. Their trend, however, is op- 
posite to ours, namely the ratio T^ LM /T™ p becomes 
smaller with increasing atomic number. This could be 
due to not taking into account the quantum character of 
the systems, which in the present paper is done via the 
factor [S + 1)/S^, where S is calculated according to 
Hund's rules. Another possible reason especially for the 
NiO result, as discussed in Ref. [lH, might be related to a 
possible importance of the short range order correlations 
that a single-site approximation like DLM would not do 
justice to. 

Concentrating on our RPA and MC results, we have to 
admit that better calculations for the individually se- 
lected TMO systems can be found in literature. Among 
them are the calculations by Zhang et a/.— for NiO (with 
a rather semi-empirical approach) and Towler et al£Q for 
MnO. However, when studying the whole TMO series 
with the same approach, such as the above DLM ap- 
plication, the work reported by Harrison^ or Wan, Yin 
and Savrasov^, it is hard to find ab initio results, treat- 
ing electron correlations at the same level of sophistica- 
tion and predicting the Neel temperatures qualitatively 
and quantitatively as accurately as in the present paper 
throughout the whole TMO series. 

To finish, we briefly discuss the Neel temperature depen- 
dence on pressure shown in Fig. |4] for all the studied 
TMOs. To calculate the pressure, p, the Murnaghan 
equation of stated has been used. Based on the be- 



haviour of J2 seen in Fig. [3J it is not surprising that 
for the whole TMO series the calculated Neel tempera- 
tures increase with pressure. Qualitatively, this agrees 
with previous experimental and theoretical results, indi- 
cating a stability of the antiferromagnetic structure up 
to high pressures (several tens of GPa, at least) before 
it collapses and a paramagnetic or low spin configura- 
tion takes over 4£J&62£M±£Z£ML™. We can compare 
the pressure dependency of T/v to experiment (for MnO 
- Ref. i, FeO - Ref. M, CoO - Ref. (H, NiO - Ref. Mj 
for p > by assuming them to be linear. Taking the 
pressure derivative of the normalized Neel temperatures, 
9(Tjv /Tm{p = 0))/dp, we find that our calculated values 
increase too slowly, roughly by a factor of 1/2. 



C. Magnon Spectra 

Considering the above results for the Jjy parameters it 
is reasonable to assume that only the nearest and next- 
nearest neighbour interactions contribute significantly to 
the magnon dispersion relation, which therefore should 
be adequately represented by Eq. (fl"3)) . For the calcula- 
tion the MFT- Ji and J2 from Tabic [Hi and the theoretical 
(calc.) magnetic moment [i from Table U were used. The 
resulting magnon spectra for all the studied TMOs, in 
the AFII structure, are shown in Fig. [5] together with the 
experimental results. Generally, the agreement between 
the calculated dispersion curves and the experimental 
observations is rather good, considering the Heiscnbcrg 
Hamiltonian used in this work — anisotropy and align- 
ment energy terms are neglected. This is also the reason 
why the calculated curves fail to reproduce the non-zero 
energies at M= (—0.5,-0.5,0.5). Besides that, minima, 
maxima and curvature arc well reproduced. Furthermore 
it can be seen that except for FeO the theoretical curves 
generally underestimate the experimental energies, which 
is due to the underestimate of the Ji parameters. The 
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x r mim' r M" 




q [2n I a] 



FIG. 5: Shown are the calculated TMO spin wave disper- 
sions together with experimental data points for MnO (black 
diamonds^, open circles^), FeO^, CoO^- and NiO^, respec- 
tively. The coordinates are cartesian and in units of 2n/a. 
The path chosen along several high symmetry lines starts at 
X= (0.25, 0.25, 0.25) and goes along [qqq] to T = (0, 0, 0), then 
along \qqq] to M= (—0.5,-0.5,0.5), and further along [00g] 
to r of the neighbouring AFII Brillouin zone, then continu- 
ing along [qqO] to M" . The inlays in the MnO panel show the 
different branches along the AFII Brillouin zone. 



relative magnitude of the peak along [qqq] varies strongly, 
as one goes through the TMO scries. This effect can be 
ascribed to the changing ratio of J^j J\. 
The qualitative agreement with previous theoretical 
works, e. g. such as that of Solovyev and Terakura^i 
is good, although not in the absolute numerical terms, 
arising from different values of the Heisenberg exchange 
parameters J^. 



V. CONCLUSION 



We have used the local self-interaction correction, im- 
plemented in the multiple scattering theory in the frame- 
work of KKR in combination with the magnetic force the- 
orem to study magnetic interactions in transition metal 
monoxides. Specifically, we have calculated the J\ and 
Ji exchange parameters, the corresponding Neel temper- 
atures and the respective magnon spectra for the whole 
TMO scries. The most important conclusion of this work 
is that the combined approach used here provides an ad- 
equate description of magnetic interactions for the series 
as a whole. Without considering correlation effects the 
theoretical results in general do not agree with experi- 
mental findings. Furthermore, we have shown that our 
ab-initio approach yields upper (MFA) and lower lim- 
its (RPA, MC simulations) for the Neel temperatures for 
the whole TMO series, and the calculated magnon spec- 
tra are in good qualitative agreement with experiment 
and other theoretical calculations. 



Acknowledgements 



We would like to thank Julie Staunton for helpful dis- 
cussions and comments. This work was supported by the 
Deutsche Forschungsgesellschaft (DFG) via the SFB 762 
"Functionality of Oxidic Interfaces" . Calculations were 
performed at the John von Neumann Institute for Com- 
puting in Jiilich, Germany. Research at the Oak Ridge 
National Laboratory was sponsored by the Division of 
Materials Sciences and Engineering, Office of Basic En- 
ergy Sciences, US Department of Energy, under Contract 
DE-AC05-00OR22725 with UT-Battelle, LLC. 



11 



1 V. V. Struzhkin, H.-k. Mao, J. Hu, M. Schwoerer-Bohning, 
J. Shu, R. J. Hemley, W. Sturhahn, M. Y. Hu, E. E. Alp, 
P. Eng, et al., Phys. Rev. Lett. 87, 255501 (2001). 

2 M. Lines and E. Jones, Phys. Rev. 139, A1313 (1965). 

3 J. R. Patterson, C. M. Aracne, D. D. Jackson, V. Malba, 
S. T. Weir, P. A. Baker, and Y. K. Vohra, Phys. Rev. B 
69, 220101 (R) (2004). 

4 C. S. Yoo, B. Maddox, J. H. P. Klepeis, V. Iota, W. Evans, 
A. McMahan, M. Y. Hu, P. Chow, M. Somayazulu, 
D. Hausermann, et al., Physical Review Letters 94, 115502 
(2005). 

5 A. Mattila, J.-P. Rueff, J. Badro, G. Vanko, and A. Shukla, 
Physical Review Letters 98, 196404 (2007). 

6 D. Kasinathan, J. Kunes, K. Koepernik, C. V. Diaconu, 
R. L. Martin, I. D. Prodan, G. E. Scuseria, N. Spaldin, 
L. Petit, T. C. Schulthess, et al., Phys. Rev. B 74, 195110 
(pages 12) (2006). 

7 U. Wdowik and D. Legut, J. Phys. Chem. Sol. 69, 1698 
(2008). 

8 A. L. Goodwin, M. G. Tucker, E. R. Cope, M. T. Dove, 
and D. A. Keen, Phys. Rev. B 72, 214304 (2005). 

9 A. L. Goodwin, M. G. Tucker, M. T. Dove, and D. A. Keen, 
Physical Review Letters 96, 047209 (pages 4) (2006). 

10 A. L. Goodwin, M. T. Dove, M. G. Tucker, and D. A. 
Keen, Phys. Rev. B 75, 075423 (pages 9) (2007). 

11 F. Ott, J. Phys.: Condens. Matter 20, 264009 (2008). 

12 P. J. van der Zaag, Y. Ijiri, J. A. Borchers, L. F. Feiner, 
R. M. Wolf, J. M. Gaines, R. W. Erwin, and M. A. Ver- 
heijen, Phys. Rev. Lett. 84, 6102 (2000). 

13 O. Bengone, M. Alouani, P. Blochl, and J. Hugel, Phys. 
Rev. B 62, 16392 (2000). 

14 A. Rohrbach, J. Hafner, and G. Kresse, Phys. Rev. B 69, 
075413 (2004). 

15 W.-B. Zhang, Y.-L. Hu, K.-L. Han, and B.-Y. Tang, Phys. 
Rev. B 74, 054421 (2006). 

16 W. M. Temmerman, A. Svane, Z. Szotek, and H. Winter, 
in Electronic Density Functional Theory: Recent Progress 
and New Directions, edited by J. F. Dobson, G. Vignale, 
and M. P. Das (Plenum, New York, 1998), p. 327. 

17 W. Temmerman, A. Svane, Z. Szotek, H. Winter, and 
S. Beiden, in Electronic Structure and Physical Properties 
of Solids - The use of the LMTO Method (Springer, Berlin 
Heidelberg New York, 2000), Lecture notes in Physics. 

18 A. Svane and O. Gunnarsson, Physical Review Letters 65, 
1148 (1990). 

19 Z. Szotek, W. M. Temmerman, and H. Winter, Phys. Rev. 
B 47, 4029 (1993). 

20 D. Kodderitzsch, W. Hergert, W. M. Temmerman, 
Z. Szotek, A. Ernst, and H. Winter, Phys. Rev. B 66, 
064434 (2002). 

21 M. Dane, M. Liiders, A. Ernst, D. Kodderitzsch, W. Tem- 
merman, Z. Szotek, and W. Hergert, Journal of Physics: 
Condensed Matter 21, 045604 (2009). 

22 C. Franchini, V. Bayer, R. Podloucky, J. Paier, and 
G. Kresse, Phys. Rev. B 72, 045132 (pages 6) (2005). 

23 X. Feng, Phys. Rev. B 69, 155107 (2004). 

24 J. Kunes, V. I. Anisimov, S. L. Skornyakov, A. V. Lukoy- 
anov, and D. Vollhardt, Physical Review Letters 99, 
156404 (2007). 

25 J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

26 A. Liechtenstein, M. Katsnelson, V. Antropov, and 



V. Gubanov, J. of Mag. Mag. Mat. 67, 65 (1987). 

27 M. Liiders, A. Ernst, M. Dane, Z. Szotek, A. Svane, 
D. Kodderitzsch, W. Hergert, B. L. Gyorffy, and W. M. 
Temmerman, Phys. Rev. B 71, 205109 (2005). 

28 I. D. Hughes, M. Dane, A. Ernst, W. Hergert, M. Liiders, 
J. Poulter, J. B. Staunton, A. Svane, Z. Szotek, and W. M. 
Temmerman, Nature 446, 650 (2007). 

29 I. Hughes, M. Dane, A. Ernst, W. Hergert, M. Liiders, J. B. 
Staunton, Z. Szotek, and W. Temmerman, New Journal of 
Physics 10, 063010 (2008). 

30 The Heisenberg Hamiltonian can be defined in several 
ways. Often the sum is multiplied with the factor 1/2, 
which corresponds to counting each y'-pair only once. 
Sometimes the minus sign is omitted. In our case also the 
absolute values of the spin vectors Si and Sj are included 
in the Jij and instead the unit vectors and are used. 
One has to take care of this when comparing exchange pa- 
rameters Jij of different works. 

31 J. Rusz, I. Turek, and M. Divis, Phys. Rev. B 71, 174408 
(2005). 

32 E. §a§ioglu, L. M. Sandratskii, and P. Bruno, Phys. Rev. 
B 70, 024427 (2004). 

33 P. W. Anderson, Theory of Magnetic Exchange Interac- 
tions: Exchange in Insulators and Semiconductors, vol. 14 
of Solid State Physics (Academic Press, New York, 1963). 

34 D. Landau and K. Binder, A Guide to Monte Carlo Simu- 
lations in Statistical Physics (Cambridge University Press, 
2000). 

35 N. Metropolis, A. Rosenbluth, M. Rosenbluth, A. Teller, 
and E. Teller, J. Chem. Phys. 21, 1087 (1953). 

36 W. Xue, G. S. Grest, M. H. Cohen, S. K. Sinha, and 
C. Soukoulis, Phys. Rev. B 38, 6868 (1988). 

37 K. Binder, Z. Phys. E 43, 361 (1981). 

38 X. Wan, Q. Yin, and S. Y. Savrasov, Phys. Rev. Lett. 97, 
266403 (2006). 

39 W. A. Harrison, Phys. Rev. B 76, 054417 (2007). 

40 For the MC simulations this factor for the Hamiltonian 
corresponds to a scaling of the temperature with the same 
factor. A quantum mechanical calculation for the MFA 



yields^ Q AB 



3| 5 -^+iij AS (0) for the TMO. So going 
from classical to quantum treatment also corresponds to 
multiplying with (S + 1)/S. Therefore, we multiplied the 
classical RPA result with (S + l)/S, analogous to the MFA. 
The values obtained from this show excellent agreement 
with such we get by using the RPA approach by Lines^, 
which is a quantum approach that, however, considers only 
nearest and next-nearest neighbour interaction and is valid 
only for the TMO and materials with the same magnetic 
structure. 

41 I. V. Solovyev and K. Terakura, Phys. Rev. B 58, 15496 
(1998). 

42 G. I. Landoldt-Bornstein, New Series, Numerical Data and 
Functional Relations in Science and Technology, vol. 27g, 
Various Other Oxides (Springer Verlag, 1992). 

43 W. Jauch and M. Reehuis, Phys. Rev. B 67, 184420 (2003). 
4 H. Fjellvag, F. Gronvold, S. Stolen, and B. Hauback, Jour- 
nal of Solid State Chemistry 124, 52 (1996). 

45 W. L. Roth, Phys. Rev. 110, 1333 (1958). 

46 W. Jauch and M. Reehuis, Phys. Rev. B 65, 125111 (2002). 

47 A. K. Cheetham and D. A. O. Hope, Phys. Rev. B 27, 
6964 (1983). 



12 



Z. Fang, I. V. Solovyev, H. Sawada, and K. Terakura, Phys. 
Rev. B 59, 762 (1999). 

49 G. Pepy, J. Phys. Chem. Sol. 35 (1974), 47. 

50 G. E. Kugel, B. Hennion, and C. Carabatos, Phys. Rev. B 
18, 1317 (1978). 

51 K. Tomiyasu and S. Itoh, J. Phys. Soc. Jpn. 75, 084708 
(2006), 43. 

52 M. D. Rechtin and B. L. Averbach, Phys. Rev. B 6, 4294 
(1972). 

53 R. Shanker and R. A. Singh, Phys. Rev. B 7, 5000 (1973). 

54 M. Hutchings and E. Samuelson, Phys. Rev. B 6, 3447 
(1972). 

55 I. P. R. Moreira, F. Illas, and R. L. Martin, Phys. Rev. B 
65, 155102 (2002). 

56 S. Shallcross, A. E. Kissavos, V. Meded, and A. V. Ruban, 
Phys. Rev. B 72, 104437 (2005). 

57 J. B. Goodenough, Magnetism and the Chemical Bond (In- 
terscience, New York, 1963). 

58 M. S. Seehra and T. M. Giebultowicz, Phys. Rev. B 38, 
11898 (1988). 

59 J. Rusz, I. Turek, and M. Divis, Phys. Rev. B. 71, 174408 
(2005). 



M. D. Towler, N. L. Allan, N. M. Harrison, V. R. Saunders, 
W. C. Mackrodt, and E. Apra, Phys. Rev. B 50, 5041 
(1994). 

F. D. Murnaghan, Proc Natl Acad Sci U S A 30, 244 
(1944). 

Y. Ding, Y. Ren, P. Chow, J. Zhang, S. C. Vogel, B. Win- 
kler, J. Xu, Y. Zhao, and H. K. Mao, Phys. Rev. B 74, 
144101 (pages 4) (2006). 

W. B. Holzapfel and H. G. Drickamer, Phys. Rev. 184, 323 
(1969). 

D. Bloch, F. Chaisse, and R. Pauthenet, Journal of Applied 
Physics 37, 1401 (1966). 

V. A. Sidorov, Applied Physics Letters 72, 2174 (1998). 
T. Okamoto, H. Fujii, Y. Hidaka, and E. Tatsumoto, Jour- 
nal of the Physical Society of Japan 23, 1174 (1967). 
D. G. Isaak, R. E. Cohen, M. J. Mehl, and D. J. Singh, 
Phys. Rev. B 47, 7720 (1993). 

J. Badro, V. V. Struzhkin, J. Shu, R. J. Hemley, H.-k. 
Mao, C.-c. Kao, J.-P. Rueff, and G. Shen, Physical Review 
Letters 83, 4101 (1999). 
M. Lines, Phys. Rev. 135, A1336 (1964). 



